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of muonic helium atom (fie^He) on the basis of a perturbation theory. We obtain 
total result for the muonically excited state hyperfine splitting Av h ^ s = 4295.66 MHz 
Q , which improves previous calculations due to the account of new corrections and more 

'_ accurate treatment of the electron vertex contribution. 

> 
OO 

■ I. INTRODUCTION 

OO 

CN| ! The muonic helium atom (fxe^He) represents the simple three-body atomic system. The 

interaction between magnetic moments of the muon and electron leads to the hyperfine struc- 
ture (HFS) of the S-wave energy levels. Whereas the ground state hyperfine splitting was 
^ ■ measured many years ago with sufficiently high accuracy in HQ Ai// p s =4465.004(29) MHz, 

the excited states in muonic helium atom were studied significantly smaller. Nevertheless, 
during the formation of muonic helium atoms the negative muon is kept by positive a- 
particle and is found to be in excited state. The lifetime of excited state with the muon 
locating in the 25-state is sufficient to measure its hyperfine splitting. This can be important 
for the check of quantum electrodynamics for three-particle muonic bound states. 

At present, hydrogen-like bound-state physics in QED is reliably developed 0-0|- Con- 
trary to two-particle bound states whose energy levels were accurately calculated in quantum 
electrodynamics j3|-0| , the hyperfine splitting of three-particle muonic states was studied on 
the basis of perturbation theory (PT) and variational approach with smaller accuracy 0- 
12]. There are several calculations devoted to HFS of the excited state ls^^s^ m muonic 
helium atom 13l-ll6|. In these papers different order corrections to HFS were studied by 
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perturbation theory, variational method and Born-Oppenheimer approach. The accuracy of 
the calculations lies in the interval of several tenth parts of MHz. The variational method 
gives high numerical accuracy of the calculation what was demonstrated in 12. 13, 15 



The correlated wave functions of the state lsfj 2 2s^ 2 were constructed and higher order 
corrections to the hyperfme splitting were also accounted for in this approach in (l5j . It is 
appropriate to mention here that three-particle kinematics was chosen for the description 
of muonic helium in @, B EMEI by different ways. So, numerical results of obtained cor- 



rections of different order in a and the ratio of the particle masses are difficult for direct 
comparison. 

Bound particles in muonic helium have different masses m e <C m M <C m a . As a result the 
muon and a-particle compose the pseudonucleus (fi 2 He) + and muonic helium atom looks as 
a two particle system in first approximation. The existence of such mass hierarchy enables 
to formulate a perturbation theory for the calculation of the energy levels. The perturbation 
theory approach to the investigation of the hyperfme structure of muonic helium based on 
nonrelativistic Schrodinger equation was developed in Using their method we describe 
the excited state in three-particle bound system (/ie 2 He) by the Hamiltonian: 



H = H + AH + AH rec , H = — Vl - — Vl - — - — , (1) 

2M ll 2M e e Xfl s fi ' 



AH = — --, AH rec =-— V -V e , (2) 



l ' >' ! i r ^ ^ 

fie -°e IIL c 



where x M and x e are the coordinates of the muon and electron relative to the helium nucleus, 
M e = m e m a / '{m e + m a ), M M = m^m a / {m^ + m a ) are the reduced masses of subsystems 
(e 2 He) + and (fi 2 He) + . The hyperfme part of the Hamiltonian contributing to HFS of S- 
states is 

3m e m /J 4 

where cr e and cr are the spin matrices of the electron and muon. In initial approximation 
the wave function of the excited state ^Sy 2 2sy 2 has the form: 

*o(xe,x,) = Mx e )M*») = ^^ W ^' 2 ( l ~ ^p) e- w ^e~ w ^% (4) 

where = 2aM M , W e = aM e . Then basic contribution to the singlet-triplet hyperfme 
splitting can be calculated analytically from contact interaction flS]): 

hfs 8tto(1+^) u f (1 + <) ( M M 2 e \ 8a 4 M e 3 
Al V =< -^T— - x e) >= —5 I 1 — TT~ + 777 ) ' = TTTTTT- ( 5 ) 



\ M J 

Numerically, the Fermi splitting is equal Vp = 4516.915 MHz. We express further the hy- 
perfme splitting contributions in the frequency unit using the relation AE h ^ s = 2nhAu h ^ s . 
For the sake of definiteness we write all numerical values with an accuracy 0.001 MHz. 
The recoil correction of first order in the ratio Mg/M^ in §5§ amounts to Au^ = 
-6j^u F = -134.769 MHz. The muon anomalous magnetic moment correction is equal 
to ^^=5.266 MHz. Modern numerical values of fundamental physical constants are taken 
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from the paper (l7| : the electron mass m e = 0.510998928(11) • 10 3 GeV, the muon mass 
m M = 0.1056583715(35) GeV, the fine structure constant a" 1 = 137.035999074(44), the 
helium mass m a = 3.727379240(82) GeV, the electron anomalous magnetic moment n e = 
1.15965218076(27) -10^ 3 , the muon anomalous magnetic moment k m = 1.16592091(63) -lO" 3 . 

In this work we aim to extend the approach of Lakdawala and Mohr [8J to the case of HFS 
of the excited state ^s^ 2 2sy 2 m muonic helium atom. The terms of the Hamiltonian AH 
and AH rec give important recoil corrections in second order perturbation theory. A feature 
that distinguishes light muonic atoms among other simplest atoms is that the structure 
of their energy levels depends strongly on the vacuum polarization (VP) effects, nuclear 



structure and recoil corrections [18H21I]. So, we investigate the role of all such potentially 
important corrections to excited state HFS in muonic helium. Another purpose of our study 
consists in the improved evaluation of the electron one-loop vertex corrections to HFS of 
order a 5 using analytical expression of the electron electromagnetic form factors 22| . 



II. RECOIL CORRECTIONS 



Along with small parameter a there are two other small parameters in this task equal to 
the ratio of the particle masses M e /M^ and M e /m a . They determine the recoil corrections 
which can be calculated by the perturbation theory. First recoil correction is presented 
in ([5]). Let us consider the recoil contribution in second order perturbation theory. For its 
study we use the initial expression for second order correction to the HFS from 0, 22| : 



Ai^'W = 16 7^ M) / &i / ^ 2 f dx 3 ^(x 3 )C (x 3 ) x (6) 

^'^,n(x3)^e,m(x 3 )^ i „(x 2 )^* m (xi) 



n,m 



Eu + E e — Eun — E e 



where the superscript (2) designates the correction in second order perturbation theory. The 
prime near the sum symbol means that the term with n = 2, m = is excluded. AH(x.i, x 2 ) 
is determined by (j2J), E^ is the Coulomb muon energy in the 2S'-state, E e is the Coulomb 
electron energy in the lS'-state. It is useful to extract in ([6]) several contributions. If the 
muon is in the 2S intermediate state (E^ — E^ n= 2 = 0) then 

Az4 /S ' (2) = 167 ? (1 + K " ] [ dx! f dx 3 |^(x 3 )|V e (x3)Ay M (x 1 )^(x 1 )G e (x 1 ,X3), (7) 
oTn e m^ J J 

where the Coulomb reduced Green's function of the electron (see (23| ) 

1 



r U v ^ - V ^."( x 3)^,n(xi) _ aMl We(xi+X3) 



,, — 'e ' J e.n 



2W eX> 



(8) 



7 \ — e '^w s x < 

ln(2H/ e x>) - ln(2H/ e x<) + Ei(2W e x<) + --2C- W e {x x + x 3 ) + 



2 ev 1 OJ 2W P x 
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x< = min(xi,Xs), x > = max{x\,x^), C = 0.577216 ... is the Euler's constant and Ei(x) is 
the exponential-integral function, 



a 



A^( Xl ) = / dx 2 C(x 2 )Aff( Xl ,x 2 )^(x 2 ) = e~ w ^ 



X\ 



1— 2„2 , l T T/3_3 



1 + 7 ™>1 + -W^ + -W^ 



(9) 

Substituting fl9]) and (JSJ) in <^ we perform analytical integration over all coordinates and 
expand exact result in the ratio M e /M^ preserving the terms of third order: 



Av. 



hfs,(2) 
2S 



+ 



Vp 1 + Ku) 



631 M e 
256 ~M U 



-14 In 



M P 5291 



M M 512 



14 In 2 



(10) 



Me 

M„ 



1391n^-^ + 1391n2 
M M 512 



61.127 MHz. 



Another correction in is determined by the muon in the IP intermediate state: 



Av. 



hfs,(2) 
2P 



167ra 2 M e 2 (l + /€ M )^ 
9m e m^ 



/ x 2 dx 2 / x 3 dx 3 ^(x 3 )V>*(x 3 ) ^ Yi 

m(-^3) ^lm 

(n 2 ) 

m=— 1 



(11) 



y rf Xl e-^ 2 / 2 e-^ 3 / 2 AiJ(x 2 , Xl )^(x 2 )^ e ( Xl ) £ ^'(111)^(113)511(^1,^3), 

m'=— 1 

where the unit vectors rij = r«/r (i = 1, 2, 3), the partial Green's function (see [24j |) 



#1 



1 1\ /e 



11^4 ' 2W / " e x > ' 2J\4W*x% m}x\ 2W eX< 2 



1 1 



(12) 

Integrating analytically in pip , we can present the expansion of final result in M e /M M up 
to terms of third order in the form: 



Aug 8 ' (2) = v F {l + k a 



837 M e 12141 M e 2 43335 M e 3 



256 M M 512 M 2 512 M 3 



-70.919 Mif*. (13) 



In addition, there exists in ([6]) the contribution corresponding to excited states of the muon: 

167T«(1 + K. 



Av hfs ' {2) 



3m e ?7V 



J dxt J dx 2 J dx 3 V>*(x 3 )^*(x 3 ) x (14) 



o 



n^2 



X 2 - Xi 



r?/> M (x 2 )V>e(Xl 



It can be calculated analytically in the leading order in Me/M^ if we replace the exact 
Green's function of the electron G e (x 3 , x^ , E^ + E e — -E^n) by free Green's function 

M e~ fe|x3 ~ Xl / 

G e ,o(x 3 , xi, E^ + E e - E^ n ) = — ^ f , b = J2M e (E^ n — E^ — E e ). (15) 

Z7T X 3 — Xi v 
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Replacing also the electron wave function by its value at the origin if) e (0) and integrating 
over the coordinate Xi in (|T4|) we obtain (see j8l. [22|): 



g-6|x 3 — xj| 

h = I dxi-j n : = 47T 

X 3 - Xi X 2 - Xi 



1 1, , 6, ,2 

---X2-X3 +-X2-X3 +... 

2 6 



(16) 



where the expansion in b/W^ ~ ^M e jM il is performed. The term 1/6 does not contribute. 
Taking second term in ([TBI) , we use further the completeness condition: 

1 

VV,n(x 3 )^,n( X 2) = 5 ( X3 ~ * 2 ) ~~ ^,2s( X 3)^,25( X 2) ~ ^ ^M,2P,m(x 3 )^,2P,m( X 2) , 
n=£2P,2S m=—l 

(17) 

where we write explicitely subscripts 25* and IP corresponding to n = 2 muonic states. 
Then the contribution of each component from (TTTT) can be calcuated analytically together 
with the term ||x 2 — x 3 |. We obtain following results corresponding to the 2S and 2P states 
in (EI): 

a hfs,(2) ,M e 2167 a hfs,(2) , M e 837 

A,/ = -„(1 + «m)]^-256". *V = Mr + /c,)^— . (18) 

Note that the order 0(M e /M^) corrections in f lT3|) and ffTB"]) corresponding to 2P-state cancel 
in the sum. To increase the accuracy of the calculation we analyze third term in ffTB"]) . It 
can be expressed in terms of the following sum jl, [s| : 

Si/2 = £ ( ^"~M 1/2 1 < ft25|^,n> | 2 . (19) 

n^2 V / ^ 

The discrete and continuum states contributions to (fT9~|) are equal correspondingly [25| : 

, v^2 16 n 6 (n 2 -l)(n-2) 2 "- 6 +i . x 



n>2 

,2 



n + 2 



cc _ f°° ? 5( k 2 , ^1/2 ( fc + 1 ) fc -fqrctj(2fc) _ i e fio 

5 1/2 -y o 2(fc + 4 ) (F+ , )6l _ e _ ¥ e - 1.563.... 

As a result the contribution of third term from the expansion ffTB"]) is equal to 

Az/ 3 ft/s ' (2) = ^(1 + (^J 5i /2 = 6.853 MHz. (21) 

Another set of recoil corrections in second order PT is determined by the perturbation 
AH rec = — ^-V^Ve. In this case the basic contribution to HFS is the following 



y- ^,n(x 3 )^, m (x 3 )^ n ^ / _ W> 2 



v nin 2 ) 
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where we substitute 



2 (2-^) (23) 



The muon 2S'-state does not contribute in fl22|) . In the case of the 2P muon intermediate 
state when i/ , ftn =^ (ii 2p, we can integrate at first over angle variables. Then we obtain the 
radial integral J °° x\dx2e~ X2 {2 — a; 2 /2)=0. So, we should consider the muon intermediate 
states with n 7^ 2 in fT22l) . Acting as in (fT4"l) we obtain: 



ht m, , 2a 3 M 2 M,W 3 W 3 



3m a m e m^n 2 



j dXij d* 2 J d X3 (l-^i) e -^3/2 e -^ 2 /2 x 



W M x 2 



;nin 2 ) ^^,n(x 3 )^,„(x2) 



n^2 



3 — 6|x3— xi 



X 3 - Xi 



Then we can integrate over xi and expand the result over small parameter b: 



(24) 



-6|x;l— Xi 



x 3 - Xi 



27r(n 2 n 3 ) 



Ax 3 x\ 2bx\ 



36 



+ 



15 



(25) 



Taking first term in square brackets we obtain two radial integrals which determine the 
contribution ({24"]) for the transition 2S — > nP (n > 2): 



■h 



,2x 
n 



e- x/2 e- x/n —Ll_ 2 (^Jx 2 dx, J 2 = R 2S {x)R nl {x)x 3 dx. (26) 



Extracting two contributions to the product of J\ and J 2 corresponding to descrete and 
continuum states we present the correction (T24"j) in the form 25| : 



Az^«(n ? 2) = u F (l + K^J^nn = 0-189 MHz, 



(27) 



_^2 l5 {n-2) 2n - 5 - 1 2n 6 {n 2 - 1) / x 



n>2 



(n + 2) 



2n+5+- 



2 15 fc(A; 2 + l)difc 



(4fc 2 + l) n /2 1-e 



2tt 
A- 



e -|arctan(2fc) = 3 2 41Q1 



Second term (— 7r(n 2 n 3 )x3) in the expansion 025]) is not dependent on n. So, we can present 
its contribution using the completeness condition: 



, „ ... Aa 3 M 2 M.W 3 W 3 



/ 2 j / i W>3 
ctx 2 / 1,(1x3 1 - 



x (28) 



W> 2 



|n 2 n 3 ) 



<5(x 2 - X 3 ) - ^,2S(X 3 )^ I2S (X 2 ) - ^,2P,m(x 3 )^,2P,m( X 2) 



m=— 1 



Only the 5-term in ( 128]) gives nonzero correction of second order in the ratio of particle 
masses: 

24 M 2 

& u W)( n ^ 2) = u F f- = 0.074 MHz. (29) 
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a b 
FIG. 1: Two photon exchange amplitudes in the electron-muon hyperfine interaction. 

Special attention has to be given to the recoil corrections connected with two-photon 
exchange diagrams shown in Fig {T] in the case of the electron-muon hyperfine interaction. For 
the singlet-triplet splitting the leading order recoil contribution to the interaction operator 
between the muon and electron is determined as follows BE, [26|: 



AKe{,; e ,2 7 (x^) = 8 2 a 2 In 3<Kx Me ). (30) 



m,, — mt m 



Averaging the potential AV^J C S ^ over the wave functions (J4]) we obtain the leading order 
recoil correction to the hyperfine splitting: 



A^e,2 7 = v F ? a ™ em » In V± = 0.812 MHz. (31) 

There exist also the two-photon interactions between the bound particles of muonic helium 
atom when one hyperfine photon transfers the interaction from the electron to muon and 
another Coulomb photon from the electron to the nucleus (or from the muon to the nucleus). 
Supposing that these amplitudes give smaller contribution to the hyperfine splitting we 
include them in the theoretical error. 



III. ELECTRON VERTEX CORRECTIONS 

In the initial approximation the potential of the hyperfine splitting is determined by fl3]). 
It leads to the energy splitting of order a 4 . In QED perturbation theory there is the electron 
vertex correction to the potential ([3]) which is defined by the diagram in Fig. [2]^a). In 
momentum representation the corresponding operator of hyperfine interaction has the form 



22 



AVi'Lt* 2 ) = - 8 "1 1 + K " ) ( ^ ) GW(tf) - 1 , (32) 



3m e m ll 

where G^(k 2 ) is the electron magnetic form factor. We extracted for the convenience the 
factor a/n from G^)(k 2 ) — 1 . Usually used approximation for the electron magnetic form 

factor G ( S{k 2 ) w G$(0) = 1 + K e is not quite correct in this task. Indeed, characteristic 
momentum of the exchanged photon is k ~ aM^. It is impossible to neglect it in the 
magnetic form factor as compared with the electron mass m e . So, we should use exact one- 
loop expression for the magnetic form factor 2?]| . Trying to improve the previous estimation 
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of the correction due to the electron anomalous magnetic moment we will use further exact 
one- loop expression for the Pauli form factor g(k 2 ) known from the QED calculation (see 
13) setting G^(k 2 )-l^g(k 2 ). 



G 



FIG. 2: The electron vertex corrections. The dashed line represents the Coulomb photon. The 
wave line represents the hyperfine part of the Breit potential. G is the reduced Coulomb Green's 
function. 



Using the Fourier transform of the potential (|32|) and averaging the obtained expression 
over wave functions (jl]) we represent the electron vertex correction to the hyperfine splitting 



as follows 22 



Au„ 



hfs 
vert 



a(l + n^)M e f m, 



2vr 2 M„ 



k 2 g(k 2 )dk(l 



Oh, ^ 



k 2 



Ml 
M 2 



+ k 2 



>>'? 

II ,7 



5.078 MHz. 



(33) 

Let us remark that the contribution (I33I) is of order a 5 . Numerical value (133]) is obtained 
after numerical integration with the exact one-loop expression of the Pauli electron form 
factor g(k 2 ) [27J]. If we use the value g(k 2 = 0) then the electron vertex correction is equal 
5.246 MHz. So, using the exact expression of the electron form factor g(k 2 ) in the one-loop 
approximation we observe small decrease of the vertex correction to the hyperfine splitting 
from I7 interaction. Taking the expression (I32p as an additional perturbation potential 
we have to calculate its contribution to HFS in second order perturbation theory (see the 
diagram in Fig. 12(b)). In this case the dashed line represents the Coulomb Hamiltonian 
AH (J2J). Following the method of the calculation formulated in previous section we divide 
again total contribution from the amplitude in Fig. [2](b) into several parts which correspond 
to the muon ground state (n = 2) and muon excited intermediate states (n / 2). In this 
way first contribution with n = 2S takes the form: 



A hfs 



W(l + re 



37r 2 m e m M 



^ f kg(k 2 )dk [ d Xl / dx 3 ^ e (x 3 )Al/ 1 (A;,X3)G e (x 1 ,X3)A^ t (xi)^(x 1 ; 

u JO J J 

(34) 



where AV r M (x 1 ) is defined by OH]) and 
AVi(/c,x 3 )= / dx 4 ^ At (x 4 ) 



/1 ok 2 1 o fc 4 

sin(fc|x 3 — x 4 |) sm(kxs) v 1 — ^w?, 



x 3 - x 4 



-Vv( x 4) 



H',f 



X3 



1 + 



fc 2 
w7, 



1 1 



(35) 
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Substituting the electron Green's function (|Sj) in ( 134)) we perform all coordinate integrations 
analytically and write (134"]) in integral form (ai = km e /W^)\ 

hfs + ( m e \ 2 (M e \ [°° kg{k 2 )a\{l-?,a 2 k 2 + 2a\k i )dk 



0.067 MHz, 



2vr 2 V«Mj V^/io (l + a 2 )*(a 2 + M 2 /M 2 )3 

a x (2a? + 8af + 8a? + 14) + — f— 2aT - 4a? + 14a a + 14(a? - l) 4 arctan —) 

M,, V a i / J 

where we expand also the integral function in the numerator in small parameter Mg/M^ and 
carry out last integration numerically. 

Second part of the vertex contribution (FigfS^b)) with n = 2P has the following general 
form: 

^ h l S t,2P = 8 f i 1 + ^ rkg(k 2 )dk [ d Xl / dx 2 / dx 3 / dx 4 ^(x 4 )^(x 3 )x (37) 

ST I ( \l* ( N Sillfc l X 4 - X s| ^e,n(x 3 )^ e * n (Xl) 

2^ ^,2P,m(x 4 )^,2P,m( X 2) i i 2^ F _ F A #( x l, X 2 )^(x 2 )^ e (Xl). 



m=— 1 n^O ' 

Then we make two subsequent integrations over coordinates x 2 and x 4 : 

ffcM = / (l - g^J = (38) 

= T^ T^T^ M + + 12 *i + 24*! - 24e^ + 24] . 

W° X\ L i i i J 



^(x 3 ) = / dx 4 x A *r w " Xi ( 1 - ) (n 1 n 4 ) "7"'^ 7 J| = (39) 



sin k 


x 4 - x 3 




x 4 - x 3 





4g7r (kx 3 coskx 3 - sinfcr 3 ) (l - 



l n l n 3j ; T4 



X3 1 + 



fc 2 
11 ' 2 

Finally, after analytical integration over coordinates xi and x 3 with the Green's function f[T2"j) 
we present the correction (1371) in the form: 

hfs 3a(l + K M )M e m2 f 00 %(£; 2 )a?(i - a 2 )d£; 



, ert) 2P 7r 2 M,(aM,) 2 7 r/ , 2 , ^^ij-iaV^ , M? 



a 2 + im + ^ K + ^§ a?+ 1 + 



(40) 

x|l2-12a^+^[75a?+215a?+209a 4 +9a?+15vrai(l+a 2 ) 4 -75(l+a 2 ) 4 aiarctanai + 156]} = 



= -0.309 Afffz. 



The expansion of integral function in the numerator of (l4"0l) is used up to terms of order 
0[M e j 'M^. In order to evaluate the vertex contribution due to excited states of the muon 
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we replace as above in (ITST) exact Green's function of the electron by free Green's function. 
Then our starting point for the calculation is related with the following expression: 



A hf s 



kg(k 2 )dk / dx 3 / dx^*^)^*^ 



sin k 


x 4 - x 3 | 




x 4 - x 3 





(41) 



61x3— xi 



n^2 



X 3 - Xi X 2 - Xi 



-^(x 2 )'0e(Xl). 



The integration over xi is done in ( 1T61) . The term proportional to 1/6 does not contribute 
because of the orthogonality of wave functions. The contribution of second term in square 
brackets in (1161) can be transformed initially by means of completeness condition: 



A hfs 
LXV vert,n±-2 



5cr 



3/2 



?>TX 2 m P m, 



71 



kg(k )dk I dx 3 / cbc 4 ^*(x 4 )?/>*(x 3 ; 



sin k 


x 4 - x 3 | 




x 4 - x 3 





/ 



c?x 2 



5(x 2 - x 3 ) - ^ )2S (x 3 )^ i2S (x 2 ) - t 



/i,2P,m 



xs)^, 



/i,2P,ml X 2j 



(42) 

|x 3 - x 2 |^(x 2 ). 



Then the integration over particle coordinates is performed analytically in all terms in (142]) . 
Remaining integration over k is carried out numerically. Omitting intermediate expressions 
for numerous integrals which are calculated in the same way as in previous corrections, we 
present final numerical result 



Au hfs 



0.992 MHz. 



(43) 



The electron vertex corrections investigated in this section have the order a 5 in the hy- 
perfine interval. Summary value of all obtained contributions in second order PT is equal 
to 0.750 MHz. Summing this number with the correction (13 3 j) we obtain the value 5.828 
MHz. It differs by a significant value 0.582 MHz from the result 5.246 MHz which was used 
previously by many; authors for the estimation of the electron anomalous magnetic moment 
contribution (see 





a 



FIG. 3: Vacuum polarization effects. The dashed line represents the Coulomb photon. The wave 
line represents the hyperflne part of the Breit potential. G is the reduced Coulomb Green's function. 
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IV. EFFECTS OF VACUUM POLARIZATION 



The vacuum effects change significantly the interaction ([I])-© between particles in 
muonic helium atom. One of the most important contributions is related with the vac- 
uum polarization (VP). The value of the electron vacuum polarization corrections to the 
hyperfine splitting is determined by the parameter equal to the ratio of the Compton wave 
length of the electron and the radius of the Bohr orbit in the subsystem (ii\He) + : M^a/m^— 
1.46727 .... It is impossible to use the expansion in a for such contributions. So, we calculate 
them performing the analytical and numerical integration over the particle coordinates and 
other parameters. The effect of the electron vacuum polarization leads to the appearance of 
a number of additional terms in the Hamiltonian operator describing (e — a)-, (fi — a)- and 
(e — /^-interactions which we present in the form jiol . I22I . 27| : 



av;j(i„) 



2a 



-2m e £,x e 



p(0 



Ve^im 2 + 1) 



(44) 
(45) 



AV^(|xe - x„|) = f- / dip{i)— e" 2 "^- (46) 



where x efJi = |x e — x M |. They give contributions in second order perturbation theory and are 
discussed below. In first order perturbation theory the contribution of the vacuum polariza- 
tion is connected with the modification of the hyperfine splitting part of the Hamiltonian 
([3]) (the diagram in Figj2](a)). It is described by the expression [HI, 19l ]: 



87ra(l + /t M ) cr 1 (T 2 a 
?>m e m^ 4 3tt J l 



2c2 



7T<5(x, 



■en) 



(47) 



Averaging the potential (14"T|) over the wave function (T4]) we obtain the following expression: 



hfs,(l) 
VP,e-/x 



q 2 (l + WjW^ 
9m e m^ 7r 2 



3u/3 roo 



-2W e X e 



tt5(x m - x e ) - 



l-X^j x e 



-11 / 1 _ W t* x ^ 1 x 



(48) 



2m e £|x AI -x ( ; 



where the superscript (1) designates the correction of first order PT and the subscript e — \i 
corresponds to the electron-muon interaction. There are two integrals over the muon and 
electron coordinates in (1481) which can be calculated analytically. As a result the correc- 
tion (1481) takes the form of one-dimensional integral in the parameter £ (a 2 = m e £/M M a): 



Au hfs,(l) 



6ttM u 1 + 



Me 

Mu 



+2al 



M e 
2 + — 
M„ 



M P 



1+ TT 1 4—^ 



M, 



M P 



M„ 



p(£R 



«2 



a 2 + 



Me 

Mu 



1 + ^ 



1 + 



M P 



M P 



M P 



M P 



- 1 



+ 

(49) 



M„ 



4 + — - 11—^ - 12 



M„ 



12 



+4a 2 



1 + 



Me 

M„ 



M 2 e Ml 
— - + 2 — - 

M 2 M 3 



M„ V M, 



- 1 



0.012 M#z. 



In the leading order this contribution has fifth order in a and first order in the electron-muon 
recoil parameter M e /M M . The similar contribution of the muon vacuum polarization to the 
hyperfine splitting is extremely small (~ 10 -6 MHz). 

Let us consider corrections of the electron vacuum polarization (j44p ~ ()46j) in the second 
order perturbation theory (the diagram in Fig. 12(b)). The contribution of the operator 
to the hyperfine splitting can be written as follows: 



Az/i 



hfs,(2) 



VP,ea 



32a 3 
9m P m, 



dxi / dx 2 / dx 3 / p(^)d^ip*(x 3 )ip* e (x 3 )x 



(50) 



n,m 



(X3 ) 1pe,m (x 3 ) % >n (x 2 ) ^* m (x X ) 



E u + E e — E, 



fi,n 



-^(x 2 )^ c (xi) 



Here the summation is carried out over the complete system of the electron and muon eigen- 
states excluding the state lsy 2 2sjf 2 . Single contribution to the sum in ([50]) is determined 
by the 2S muon intermediate state because of the wave function orthogonality condition: 



Ai/ 



hfs,(2) 
VP,ea 



32« 3 (1 + K p 

9m P m u . 



p(O d f / rfx i / ^x 3 |^ At (x 3 )| 2 V' e (x3)G' e (xi,X3)- 



Using the electron Green's function (JS 
Xi and x 3 and present (I5~TI) as follows: 



Az/ 



hfs,(2) 
VP,ea 



-v F - 



q( l + /y)M e 
37rM„, 



-^e(Xl). 

(51) 

we perform analytical integration over coordinates 
4 + a 2 (a 2 + 2) 2 



a 2 (a 2 + l) 4 



-0.028 M#z. 



(52) 



This contribution has the same order of the magnitude 0(a 5 M e /M^) as the previous 
one in first order perturbation theory. The similar calculation can be performed in the 
case of muon- nucleus vacuum polarization operator ( l45l) . The intermediate electron state 
is the IS state and the reduced Coulomb Green's function of the system transforms to the 
Green's function of the muon. The correction of the operator (|45p to the hyperfine splitting 
is obtained in the following integral form: 



Au hfs,(2) 



where 
G /1 (x 2 ,x 3 ) 



32a 3 (1 + k a 

9m e Tn^ 



p(£)d£ / dx 2 / dx 3 ^*(x3)Vv(x 2 )|^ e (x 3 )| 2 G M (x 2 ,x 3 ), 

(53) 



OlM 2 , e -(^+x 3 )/2 



16tt 



X2X3 



8a; < - Ax\ + 8x> + Ylx < x > - 2<ox\x > + 2x 3 < x > - Axl, - (54) 



26x < Xy+23xl : xl-xl : x^+2x < x 3 > -x\xl > +4:e x (l-x < )(x > -2)x > +4:(x < -2)x 



•2)x>x 



[-2C + £z(x<) - ln(x<) - ln(ac>)] 
Integrating over coordinates we obtain: 

a(l + K^)M ( r °° A - 4 1 «- 3 1 A ~ 2 



Au hfs ' {2) 



-v F - 



2vrM„, 



/ j-\ ,A<4 + 8a 3 + 4a 2 + 6a 2 + 1 
(1 + a 2 ) b 



-0.099 MHz. (55) 
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The vacuum polarization correction connected with the operator ( 146]) in the second order 
perturbation theory is the most difficult for the calculation. Indeed, in this case we have to 
consider the intermediate excited states both for the muon and electron. We have divided 
total contribution into several parts. The first part in which the intermediate muon is in 
the 2S state can be written as: 



2 poo 



p(Od£ / xfe-^dx! x 2 3 dx 3 (l - ^)V* 8 < 1+ ^ 

(56) 

Me 

M e (x l +x 3 ) 1-e 1 



x A W*i) [l^T-ln ^>-ln £ x<+E i(£ X< ) + 7 --2C 



Me 



M U X > 



M U X < 



-]■ 



where the function Vyp ^(^i) is equal to 
AVVp^(xi) = 



x\e~ X2 [ 1 - 



X 2 



dx<? / dz 



e -a 2 |xi-x 2 



Xl - X 2 



(57) 



2x 1 {a 2 -l) 4 



8e~ a2a;i (2c4 + 3a^ + 1) + e~ Xl a^i(x a - 2) 2 + a£(xi(10 - 3(xi - 2)x 1 ) - 16) + 



+a^(3x? - 8xi - 24) - xi(xi(xi + 2) + 6) - 8) 
Substituting f loT)) in f l5"Sj) we obtain the following result with an accuracy 0(M e /M^): 



^ U VP,ne,2S - ~ V F 



aM e (l + k. 



768vrM M J 1 a 2 (l + a 2 ) 



■(j7a 7 2 + 616al + 2151^ + 4272t4+ (58) 



+5267^ + 4168^ + 1929a 2 + 512 = -0.009 MHz 



The second part of the vacuum polarization correction to the hyperfine splitting due to the 
2P muon state can be presented as 



L ^ u VP,ne,2P 



16a 3 (l + k a 
9m e m M 



x 



n^O 



^ e , n (x 3 )^* (Xi 



dXi / dx 2 / dx 3 ^*(x 3 )^*(x 3 ) ^,2P,m(x 3 )^* i2 P,m( X 2)x 
^ m=-l 

(59) 



-Ep — K 



e,n J 1 



e -2m e g|x 2 -x 1 

p(0^— — r^^(x 2 )^ e (xi; 



x 2 — Xi 



Averaging over the directions of the vector n 3 < (n 3 n 2 )(n 3 ni) >=(nin 2 )/3 and evaluating 
the integral 



AVVr»(xi 



x\e~ X2 [ 1 - 



X 2 



dx?, I zdz- 



3 -a 2 |xi-X2 



1 



e~ Xl x (60) 



xi-x 2 | x\{a\ - l) 4 
2A(al + l)(x 1 + l)-12xj(4-l) + 2x 3 1 (al-3a 2 2 + 2)-x^(a 2 2 -l) 3 -24e- a2 * 1 (a^ + l)(a 2 x 1 + l)}, 



we find this contribution numerically: 



Al hfs,(2) _ 



9a(l + /v)M e [°° p(Z)d£ 



1927rM„ 



[i + a 2 y 



bal + 40a^ + 139a* + 272^+ (61) 
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+3230 2 . + 232a 2 + 93 = -0.019 MHz. 
Finally, there exists the contribution of the muon intermediate states with n ^ 2S, 2P: 

A^ 2 gp = - 16 Q 3 j t 1 e ^ ) / dX! J ^X 2 j rfx3^(x 3 )^(x3)g^n(x3)^, n (x 2 )x (62) 

roc e -2m e g|x 2 -xi| 

G , e (xi,X 3 ,E / , + E e - E^n) / p(£)d£— j— ^(X 2 )^ e (xi). 

J\ | x 2 — Xi| 

Then we have replaced in (|62|) exact electron Coulomb Green's function G e by free electron 
Green's function ffl5l) which contains the parameter 6 = [2M e (E^ M — E^ — E e ] l l 2 . We also 
replace the electron wave functions by their values at the origin neglecting higher order recoil 
corrections. After that the integration over xi can be done analytically: 

e -6|x3-xi| e -2m e §|x 2 -xi| ^ n e -6i|x 3 -x 2 ] _ e -a2|x 3 -x 2 

rfx i1 f - i f~ = ~7F7 1 VTTa 2^ — = ( 63 ) 

|x 3 -xi| |x 2 -xi| Wf, |x 3 - x 2 |(6f - a 2 ) 



Air 



1 - e - a2 l X3 - X2 l 2&i 

(22|x 3 — X 2 J & 2 



where we have performed small-parameter b\ = b/W^ expansion. For the further trans- 
formation we use the completeness relation (fTTl) dividing total correction (162!) into three 
terms. The coordinate integration in each of the terms can be presented individually in 
analytical form. As a result the correction (162 j) is reduced to one-dimensional integral over 
the parameter £ which is evaluated with high accuracy numerically: 

a hfs,{2) ar(l + « M )M B f°° f nKa 6 . 17Q1 5 . 4, (tiA \ 

AU ^ P = VF 384irM„ / (TT^F I 256 " 2 + 173102 + 4949fl2+ (64) 



7686a^ + 6874a 2 , + 3375a 2 + 665 ) = 0.076 MHz. 



i 2 1 uui-*u, 2 



We preserve in (1641) only leading order term in the ratio M e /M M . In whole, in (1641) we 
have second order correction in two small parameters a and M e /M M . We keep in (I63p only 
first term in the square brackets because other terms give recoil corrections of higher order. 
Numerical values of obtained corrections are presented in Table I. 

There exists another contribution of the second order perturbation theory in which we 
have the vacuum polarization perturbation connected with the hyperfine splitting part of 
the Breit potential (JTJ) (see Fig. [5]). Other perturbation is determined by AH (T5]). We can 
divide this correction into three parts. One part corresponds to the muon in the 2S'-state. 
The 5-function term in fT4Tj) gives the following contribution at n = 2 (compare with fflQ]) ): 

Av!£>® = v F (l + «„) — r P&W^r-- ( 65 ) 
vp,u,2S *\ M; 3?r y i f\v t - 2 56M fl K ' 

Obviously this integral is divergent. So, we have to consider the contribution of second term 
in (I47p to the hyperfine splitting which is determined by more complicated expression: 

16a 2 (l + K M )™ 2 f°° „ f £\ t 2 At f A „ L , f „ M2 f j„ .,. ,„ N e- 2me l X3 " X4 l 



A^fe ( 2 % = - 16a Q (1 + WMK r P (0e^ [ rfx 4 |^(x 4 )| 2 / ab^ e (x 3 ) 

(66) 



x 3 - X4 



15 




FIG. 4: Vacuum polarization effects in second order perturbation theory. The dashed line repre- 
sents first part of the potential AH (3). The wave line represents the hyperfine part of the Breit 
potential. 



x 2^ p p A^ M (x 1 )V e (x ly 



where the matrix elements entering in ( 166]) are denned by ([9]) and ( 1571) . Integrating over all 
coordinates in (166]) and summing it with (165]) we obtain the following convergent integral in 
£ in the leading order with respect to the ratio [M e jM^)\ 



a fc/a,(2) , A Ms,(2) 
^"VP11,2S ~r *-*^VP,12,2S ~~ ^ F 



a(l + K v )M e 



3vrM„ 



P(0 ^ '51a| + 4084+ (67) 



i (a 2 + l) 8 

+1177^ + 1872a^ + 2029a^ + 1464a 2 + 63l) = 0.004 MHz. 



G 



FIG. 5: Nuclear structure effects in second order perturbation theory. The bold point represents 
the nuclear vertex operator. The wave line represents the hyperfine part of the Breit potential. 

Next correction arises from the muon in 2P intermediate state. Making the same division 
of total contribution into two parts as in fl67|) we present here final result: 



A hfs,(2) Au hfs,m 



VP,12,2P 



-V F 



3a (1 + K M )M e 



2567rM M J 1 (a 2 + 1) 



5a b 9 + A0a 5 9 + 139a 



(68) 



+272^ + 323a^ + 232a 2 + 93 
Let us consider now the terms with n ^ 2S, 2P: 
16a 3 (l + K u )M t 



-0.005 MHz. 



A hfs,(2) 
£ * U VP,iM-e,rv£2S,2P 



187r 2 m P m,, 



^ P(0<% j dx l j d *2 J ^ X 3 J rfX4^;(x 4 )^ e *(x 3 )x 

(69) 
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-6|X3— XI I 



X 3 - Xi X 2 - Xi 



r ^u(x 2 )^ e (xi) 



vr5(x 3 -x 4 ) 

x-_> — x 1 1 J 

n^2S,2P 

where free Green's function of the electron ( fl5l) is used. Changing i/j e (xi) —> ip e (0) and 
using first nonzero term in the expansion (1161) we present the contribution of the 5-function 
from ([69]) in the form: 



a hfs,(2) _ _ 

^ U VP,ll,n^2S,2P ~~ U F 



a(l + K„)M e 
3ttM u 



665 
128' 



(70) 



As expected the integral (jTOj) in £ is divergent. We should consider together with it another 
term from 



A hfs,(2) _ _ 

^ A ' / yP,12,n2S , ,2P ~~ "F 



2a 2 (l + Kll )M e m 2 e 
3tt 2 



2 />oo 



p{£)d£ / rfx 2 / dx 3 / rfx 4 



3 -fe|x 3 -x 4 | 



| x 3 — x 4 



■^(x 4 )x 
(71) 



<5(x 2 - X 4 ) - ^(X 2 )^*(X 4 ) - ^ ^,2Pm(x 2 )^* i2 p m (x 4 ) 



171= — 1 



|x 3 - x 2 |^(x 2 ) 



All three terms in square brackets give the contributions in (1711) . They can be evaluated in 
the same way as in previous sections. The sum of expressions (1701) and (ITTil is determined 
by finite integral which is calculated numerically: 



a hfs,(2) a h/s,(2) 

ZAZ/ Vpil,n^2S,2P ^VP^nyspP ~ V F 



q(l + /y)M e 
3vrM„ 



+403(4 + 2075aS + 4596a^ + 7310c4 + 6559(4 + 4511a 2 , + 1639a 2 + 256 ) = 0.012 MHz. 



128a 2 (a 2 + 



r_ 5a 9 + l 6a 8+ (72 ) 



All vacuum polarization contributions to the hyperfine splitting are presented in Table I. In 
their calculation we use systematically the expansion in M e /M M retaining first order terms. 
So, all corrections evaluated in this section are of second order in two small parameters a 
and M, ;M,,. 



V. NUCLEAR STRUCTURE CORRECTIONS 

In the leading order in a the nuclear structure corrections to hyperfine splitting are 
determined by the charge radius of a-particle r a . Considering the interaction between the 
muon and the nucleus we can present the nuclear structure correction in the interaction 
operator in the form 

AV str ^ a {r^) = ^nZa < r 2 a > S(r„), (73) 

where the subscript str designates the structure correction. The contribution of the operator 
AV s tr,fi-a t° the hyperfine splitting appears in second order perturbation theory (see the 
diagram in Figj5]). We can write it in the integral form: 

A hfs 647T 2 a 2 (l + K^rl Wl 12 f , |2 

A ^str,„-a = o " 7-7E= / (x 3 )|<0 e (x 3 )| G^(x 3 , 0, EJ. (74) 

ym e m„ Zn J 
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TABLE I: Hyperfine singlet-triplet splitting of the excited state in muonic helium atom. 



Contribution to the HFS 


Av h f s , MHz 


Reference 


The Fermi splitting 


4516.915 


(5) 


The muon anomalous magnetic 
moment correction a^vp 


5.266 


(5) 


Recoil corrections in M e /M^ 
in first order PT 


-132.608 


(5) 


Recoil correction in M e /M^ in second 
order PT, muon IS state 


61.127 


(10) 


Recoil correction in M e /M^ in second 
order PT, muon IP state 


-70.919 


(13) 


Recoil correction in M e /M^ in second 
order PT, muon excited states 


-109.978 


(18), (21) 


Recoil correction of order — .,/-&- a 4 in second 

orrlpv PT* iniinTi P"ypit,pn ^tatp^ 


0.189 


(27) 


Recoil correction of order M f. a 4 in second 
order PT, muon excited states 


0.074 


(29) 


Recoil correction of order 
a b (m e /m^) m(m e /m M ) 


0.812 


(31), [15] 


Recoil correction In — vp 
in third order PT 


19.434 


M 


Electron vertex correction of order a 5 
in first order PT 


5.078 


(33) 


Electron vertex correction of order a 5 
in second order PT 


0.750 


(36), (40), (43) 


One-loop VP contribution in I7 
(efi) interaction of order a 5 


0.012 


(49) 


One-loop VP contribution in the 
electron-nucleus interaction in the 
second order PT of order a 5 


0.028 


(52) 


One-loop VP contribution in the 
muon-nucleus interaction in the 
second order PT of order a 5 


0.099 


(55) 


One-loop VP contribution in the 
eiecLi on-muon mieiaciion 1x1 me 
second order PT of order a 5 


0.059 


(58),(61),(64),(67), 


Nuclear structure correction of order 
a 6 in second order PT 


-0.015 


(75),(76) 


Relativistic correction of order a 6 


-0.050 


[15] 


Electron vertex correction of order a 6 


-0.615 


[6, 28-30] 


Summary contribution 


4295.658 
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Using the muon Green's function with one zero argument we make analytical integration 
in (I74p and present final result in the form: 

, , 8a 2 M 2 r 2 (l + k u ) ( M P M 2 \ 
A^_ Q = -v F »-f ^ \6j£ - 41^ + . . .J = -0.014 MHz. (75) 

Numerical value is obtained by means of the charge radius of the a-particle r a = 1.676 fm. 
The same approach can be used in the study of the electron-nucleus interaction. The elec- 
tron feels as well the distribution of the electric charge of a particle. The corresponding 
contribution of the nuclear structure effect to the hyperfine splitting is the following: 



Aa 2 M 2 rl \\ M P M 2 ( 3 M P , 

-TT^-S + )+■■■ = -0.001 MHz. 76 

2M„ ill 2 \2 M ' 



We have included in Table I the total nuclear structure contribution which is equal to the 
sum of numerical values f!75|) and (ITS"]) . 



VI. CONCLUSION 

In this work we have performed the analytical and numerical calculation of some basic 
contributions to the hyperfine splitting of the excited state lsj^s^ i n muonic helium 
atom connected with the recoil effects, electron vertex corrections, vacuum polarization and 
nuclear structure effects. We used the perturbation theory for this task as it was formulated 
previously in As was shown by Lakdawala and Mohr that the ground state hyperfine 
splitting in muonic helium atom can be expressed analytically by a series in the ratio M e /M^. 

We extend their method to the case of the excited state ls^J 2 2sy 2 accounting for not only 
purely recoil corrections but also other vacuum effects controlled by usual small parameter 
a in QED. Our results are presented in Table I in which we give also references to the 
calculations of some other corrections which are not considered here. We have included in 



the total value of the hyperfine splitting the recoil logarithmic correction obtained in [14 



in third order PT, the relativistic correction obtained in [15[, the electron vertex correction 
of order 0(a 2 ) known from the theory of hydrogenic atoms 0, [28 30| . Our numerical 
result 4295.66 MHz is, in whole, in the agreement with the earlier performed calculations 
by the perturbation theory and variational method: 4291.50 MHz Q, 4287.01 ± 0.10 MHz 



15] . Our analytical expressions of different recoil corrections are in agreement with the 



calculation in [l4|. Comparing our results (EEU]), ([12]), (USD, (EU), fl2H]) with [14J we can find 



their coincidence in the leading order in m e jm^ m e /m a and logarithmic terms In— . 
We should take into account that the used definition of reduced masses M e and M^ of the 



particles is different from m* and m* in [14]]. The correspondence between our results and 
basic expression (62) from [l4| can be revealed if we use the relation: M e /M^ = m* e /m^ + 
m*/m a . It is interesting to note that the recoil correction of order m e /m a appears in our 
kinematics from the definitions of reduced masses M e and M^. The recoil Hamiltonian 
AH rec gives the contribution in second order perturbation theory of order M 2 /m a M^ (129]) . 
Numerical differences of our results from [l4| are determined by different choice of three- 
particle kinematics, higher order corrections M 2 /M 2 which we account for in intermediate 
expressions and calculated vacuum polarization effects of order aM e /M^. So, for example, 
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numerical value of the recoil correction of the leading order ^— vp in fl5]) is equal to 

— 134.769 MHz. This number differs by more than 2 MHz from the value —132.608 MHz 
presented in Table I because we preserve corrections of higher order presenting fl5]) in the 
form u F (l + - QM e /M^ + 21M 2 /M 2 - 55M e 3 /M^). Another important source of the 
difference is specified by more accurate calculation of the electron vertex correction in section 
III related with the Pauli form factor. We use exact expression for it making numerical 
integration in corresponding amplitudes. Finally, we use in hyperfine splitting part of the 
Hamiltonian the factor (1 + k^) related with the muon anomalous magnetic moment what 
also gives definite modifications of numerical results. Our summary value of the vacuum 
polarization effects 0.198 MHz agrees with the result 0.229 MHz obtained by variational 



method calculation in 15 



Effects of the nuclear structure give negligibly small contribution to hyperfine splitting. 
In muonic helium atom the muon strongly interacts with a-particle at distances of order 
r ~ 1/M^a. It was discovered recently in the muonic hydrogen Lamb shift experiment that 
the experimental value differs from theoretical one on 0.3 meV with CODATA value of the 



proton charge radius [31] . This disagreement between the proton charge radius extracted 
from Lamb shift measurements of muonic and electronic hydrogen can be considered as the 
appearance of new muon physics [32]. The plans to measure the (2S — 2P) transition fre- 
quencies in muonic helium ions (/4-ffe) + and (//fiie) + with an accuracy of 50 ppm can help 
to solve the proton radius puzzle [33]. In the case of muonium the most precise compari- 



son between theory and experiment is possible [6J, |34j . It will be interesting to use muonic 



helium atom for the comparison of theoretical value of hyperfine splitting 4295.66 MHz 
with experimental data for the revealing of new interactions in muonic sector. Main the- 
oretical error in our calculation is connected with uncalculated in full recoil contribution 
of second order vftA ~ 0.11 MHz. In addition, the electron vertex corrections in two- 

photon processes which are estimated in this work very approximately also give essential 
error upa 2 ~ 0.25 MHz. Thereby, the total theoretical error does not exceed ±0.40 MHz, 
but large values of coefficients can essentially change preliminary estimate. Further im- 
provements of the theoretical result require more careful construction of the three-particle 
interaction operator connected with the multiphoton exchanges. 
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